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DESCRIPTION 
MAGNETIC RESONANCE IMAGER 

5 Technical Field 

The present invention relates to a magnetic resonance 
imaging technology. 

Background Art 

10 Magnetic resonance imaging (MRI) apparatuses are 

medical-purpose image diagnosis systems that cause hydrogen 
atomic nuclei, which are contained in any transverse plane of 
a subject, to exhibit nuclear magnetic resonance, and produce 
a tomographic image of the plane using an induced nuclear 

15 magnetic resonance signal. 

In general, when a slicing magnetic field gradient to 
be used to select a subject's plane whose tomographic image 
is to be produced is applied, an excitation pulse that excites 
magnetization in the plane is applied at the same time. Thus, 

20 a nuclear magnetic resonance signal (echo) is induced in the 
stage of convergence of the excited magnetization. In efforts 
to append positional information to magnetization, a 
phase-encoding magnetic field gradient and a readout magnetic 
field gradient that are perpendicular to each other in a section 

25 are applied until an echo results from excitation is acquired. 
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Ameasured echo is mapped to a k-space having an axis of abscissas 
kx and an axis of ordinates ky. Image reconstruction is then 
performed through inverse Fourier transform. 

A pulse with which an echo is induced and each magnetic 
5 field gradient are applied based on a predefined pulse sequence . 
As for the pulse sequence, various pulse sequences intended 
for respective purposes are known. For example, a gradient 
echo (GE) -type fast imaging method is a method in which a pulse 
sequence is repeatedly applied, a phase-encoding magnetic field 

10 gradient is sequentially changed for every repetition in order 
to measure the number of echoes required for producing one 
tomographic image . 

Fig. 1(A) shows a pulse sequence to be employed in GE 
radial scanning (refer to, for example, ^^Magnetic Resonance 

15 Imaging - Physical Principles and Sequence Design" by E. Mark 
Haacke, et al. (Wiley-Liss, pp. 303-330, 1999)). Actions to 
be performed for the pulse sequence will be described below. 

Along with application of a z-direction slicing magnetic 
field gradient 201, a radiof requency (RF) magnetic field pulse 

20 202 for excitation of magnetization of protons at a resonant 
frequency f 0 is applied in order to induce a nuclear magnetic 
resonance phenomenon in protons in a certain slice of an object 
entity. After dephasing magnetic field gradient pulses 203, 
204, and 205 are applied, while readout magnetic field gradient 

25 pulses 206 and 207 are applied, a nuclear magnetic resonance 
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signal (echo) 208 is measured. After the measurement of an 
echo is completed, re-phasing magnetic field gradient pulses 
209, 210/ and 211 are applied in order to restore the phase 
of magnetization in preparation for the next excitation. 
5 The above procedure is repeated Ne times with a repetition 

time TR determined, whereby Ne echoes are measured. The 
dephasing magnetic field gradient pulses 204 and 205, readout 
magnetic field gradient pulses 206 and 207 , re-phasing magnetic 
field gradient pulses 209 and 210 have the amplitudes thereof 

10 changed for every repetition as shown in Fig. lA. In the case 
of the illustrated sequence, the dephasing magnetic field 
gradient pulse 204 and re-phasing magnetic field gradient pulse 
209 change step by step from -Ne/2 to Ne/2-1. The dephasing 
magnetic field gradient pulse 205 and re-phasing magnetic field 

15 gradient pulse 210 change step by step from 0 through -Ne/2 
to -1. The readout magnetic field gradient pulse 206 changes 
step by step from Ne/2 to -Ne/2-1. The readout magnetic field 
gradient pulse 207 changes step by step from 0 through Ne/2 
to 1. 

20 Measured echoes are, as shown in Fig. 1(B), mapped to 

a k-space . The drawing is concerned with a case where Ne denotes 
128 . In the k-space, one echo is expressed with one line passing 
through an origin O, and echoes are disposed equidistantly in 
a rotating direction. A difference AG between the angles of 

25 adjoining echoes is 7i/Ne radian. 
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The k-space is transformed into a Cartesian grid by 
performing gridding (refer to, for example, ''Selection of a 
Convolution Function for Fourier Inversion Using Gridding'' by 
Jackson JI, Meyer GH, Nishimura DG (IEEE Trans. Med. Imaging, 
5 Vol. 10, No. 3, pp. 473-478, 1991)). Thereafter, image 

reconstruction is performed through two-dimensional inverse 
Fourier transform. An imaging time required for one image 
corresponds to a product of a TR by the number of echoes. For 
example, assuming that one image is reconstructed using one 

10 hundred and twenty-eight echoes with the TR set to 4 ms, the 
imaging time comes to 512 ms. 

In order to reconstruct an image having N pixels in rows 
and columns, the number of samples per echo and the number of 
echoes are normally set to N. If the number of echoes is smaller 

15 than N, the imaging time is shortened and a temporal resolution . 
improves. For example, assuming that only odd-numbered echoes 
shown in Fig. 1(B) are measured, the number of echoes is 64 
and the imaging time is a half of the above imaging time. 

20 Disclosure of the Invention 

Problem to be Solved by the Invention 

However, when the temporal resolution is improved 
according to the foregoing method, a spatial resolution is 
degraded and artifacts occur. In this case, the k-space is 

25 expressed as shown in Fig. 2. In the drawing, a dot line 215 



signifies that an echo has not been measured. Compared with 
the k~space shown in Fig. 1 (B) , the k-space is short of samples. 
Therefore, the spatial resolution of a reconstructed image is 
degraded and artifacts occur . Namely, when the number of echoes 
5 is decreased in order to improve the temporal resolution, 
artifacts occur to deteriorate image quality. 

Moreover, the procedure starting with application of a 
slicing magnetic field gradient pulse and ending with 
measurement of an echo like the one shown in Fig. 1 (A) is repeated 

10 with a repetition time TR determined appropriately. Echoes 
required for one image are thus measured. The echoes are mapped 
to the k-space as shown in Fig. 1(B), and then, as shown in 
Fig. 8(B) , re-mapped to a grid-like k-space 802 by performing 
gridding. Thereafter, two-dimensional inverse Fourier 

15 transform is performed for image reconstruction. An imaging 
time required for one image is 0.256 sec on the assumption that 
the TR is set to 4 ms and sixty-four echoes are acquired. For 
imaging of the heart, since a cardiac cycle is approximately 
1 sec, the cardiac motion is a non-ignorable factor affecting 

20 the image quality. If imaging is not performed in real time, 
the adverse effect of the cardiac motion can be suppressed 
according to a method such as cardiac gating. As for real-time 
imaging, a subject makes a large motion during imaging, and 
artifacts occur. 

25 An object of the present invention is to provide a magnetic 
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resonance imaging technology for efficiently suppressing 
artifacts in radial scanning. 



Means for Solving the Problem 
5 In order to accomplish the aforesaid object, a magnetic 

resonance imager in accordance with the present invention has 
features described below. 

1 . The magnetic resonance imager includes a control unit 
that controls a pulse sequence according to which a 

10 radiofrequency magnetic field and magnetic field gradients are 
applied to a subject lying down in a static magnetic field in 
order to detect a magnetic resonance signal induced in the 
subject, and a processing unit that handles the signal. The 
control unit (1) controls a pulse sequence based on which radial 

15 scanning is achieved, (2) applies the pulse sequence a plurality 
of times so as to acquire image echoes, and (3) applies the 
pulse sequence a plurality of times so as to acquire reference 
echoes each of which lies among image echoes in the k-space. 
The processing unit (1) divides the image echoes and reference 

20 echoes into a plurality of groups, (2) uses a reference echo 
and preceding and succeeding image echoes to calculate an 
estimation coefficient, (3) uses the estimation coefficient 
to estimate unmeasured echoes lying among the image echoes in 
the k-space. 

25 2 . The magnetic resonance imager includes a control unit 
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that controls a pulse sequence according to which a 
radiof requency magnetic field and magnetic field gradients are 
applied to a subject lying down in a static magnetic field in 
order to detect a nuclear magnetic resonance signal induced 
5 in the subject. The control unit (1) detects the nuclear 
magnetic resonance signal by radially scanning a k-space, (2) 
produces a plurality of images, (3) adopts a sliding window, 
and (4) performs the scanning at intervals of n echoes. Thus, 
artifacts are suppressed by a temporal filter, 
10 According to the present invention, there is provided 

a magnetic resonance imager capable of efficiently suppressing 
artifacts in radial scanning. 

Brief Description of the Drawings 
15 Fig. 1 is an explanatory diagram concerning a pulse 

sequence for conventional GE radial scanning and a k-space; 

Fig. 2 is an explanatory diagram showing a conventional 
k-space employed for radial scannings- 
Fig . 3 shows an example of the configuration of a nuclear 
20 magnetic resonance imager to which the present invention is 
adapted; 

Fig . 4 is an explanatory diagram showing a k-space employed 
for radial scanning according to the present invention (first 
embodiment) ; 

25 Fig. 5 is an explanatory diagram showing arrangement of 



echoes employed for estimation of unmeasured echoes according 
to the present invention (first embodiment); 

Fig. 6 is a flowchart describing estimation of unmeasured 
echo to be performed according to the present invention (first 
5 embodiment) ; 

Fig. 7 shows the results of estimation of unmeasured echo 
performed according to the present invention (first 
embodiment) ; 

Fig. 8 is an explanatory diagram concerning a sequence 
10 of scanning a k-space at intervals of N echoes employed in the 
present invention (second embodiment); 

Fig. 9 shows frequency-components of a motion picture 
reconstructed by scanning the k-space at intervals of three 
echoes, and a temporal filer (low-pass filter) (second 
15 embodiment) ; 

Fig. 10 is an explanatory diagram concerning the results 
of imaging performed according to the present invention (second 
embodiment) ; and 

Fig. 11 is an explanatory diagram concerning an imaging 
20 procedure employed in the present invention (second 
embodiment) . 
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Best Mode for Carrying out the Invention 

Referring to the drawings, embodiments of the present 
invention will be described below. 
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(First Embodiment) 

a' 

Fig. 3 is a block diagram showing the outline configuration 
of a magnetic resonance imager. 

In Fig. 3, reference numeral 101 denotes a magnet that 
5 generates a static magnetic field, 102 denotes a coil that 
induces a magnetic field gradient, and 103 denotes a subject 
(for example, a living body) . The subject 103 is carried to 
a space of a static magnetic field generated by the magnet 101. 
A sequencer 104 transmits an instruction to each of a gradient 

10 power supply 105 and a transmitter 106, and thus allows them 
to induce a magnetic f ieldgradient and a radiof requencymagnetic 
field respectively. The radiof requency magnetic field is 
applied to the subject of examination 103 via a probe 107. A 
signal induced in the subject of examination 103 is received 

15 by the probe 107, and detected by a receiver 108. A nuclear 
magnetic resonant frequency to be adopted as a reference for 
detection (hereinafter a detection reference frequency) is set 
by the sequencer 104. The detected signal is transmitted to 
a computer 109. The computer 109 performs signal processing 

20 such as image reconstruction. 

The results of signal processing are displayed on a display 
110. If necessary, a detected signal and the conditions for 
measurement may be stored in a storage 111. Moreover, an 
electrocardiograph 114 connected to the sequencer 104 is located 

25 in the static magnetic field space, and can be used to obtain 
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an electrocardiogramof the sub j ect 103 . The electrocardiogram 
is transferred to the sequencer 104. If the homogeneity in 
a static magnetic field should be adjusted, a shim coil 112 
is employed. The shim coil 112 supports a plurality of channels, 
5 and a shim power supply 113 supplies power to the shim coil 
112. When the homogeneity in the static magnetic field is 
adjusted, a current flowing through each shim coil is controlled 
by the sequencer 104. The sequencer 104 transmits an 
instruction to the shim power supply 113, and thus allows the 

10 coil 112 to induce an additional magnetic field with which the 
inhomogeneity in the static magnetic field is corrected. 

The sequencer 104 extends control so that the other 
components will operate at respective preprogrammed timings 
to respective preprogrammed extents . Among programs , a program 

15 stating the timing of signal detection and the timings and 
amplitudes of a radiofrequency magnetic field andmagnetic field 
gradients is referred to as a pulse sequence. 

In the present embodiment, a GE sequence shown in Fig. 
1 is adopted as a pulse sequence. A TR for the pulse sequence 

20 is 4 ms, and a number of repetitions comes to 72 by adding eight 
repetitions for reference to a conventional number of 
repetitions of 64 . Thus, seventy-two echoes are measured. The 
positions of echoes in the k-space are indicated with solid 
lines in Fig. 4. 

25 In Fig. 4, among echoes indicated with solid lines, an 
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echo 213 indicated with a boldface line is a reference echo. 
The echo is used to calculate an estimation coefficient, whereby 
unmeasured echoes 212 indicated with dot lines are estimated. 
The sum total of echoes including estimated ones comes to one 
5 hundred and twenty-eight echoes. The procedure will be 
described in conjunction with the flowchart of Fig. 6. 

To begin with, one hundred and twenty-eight echoes 
including the unmeasured echoes 212 are divided into eight groups 
one of which includes seventeen echoes (step 401) . As shown 

10 in Fig . 5 , marginal echoes of each group are identical to marginal 
echoes of two adjoining groups. In Fig. 5, an echo 214 is 
included as a marginal echo in each of the first group 301 and 
second group 302. The reference echo 213 is measured to be 
included in the middle of each group. 

15 Thereafter, an estimation coefficient A=[al,a2] to be 

used to estimate the unmeasured echoes 212 included in each 
group is calculated according to the expression (1) below (step 
402) . 

A=RS- (1) 
20 where S equals [SI, S2] "T (a"T denotes a transposed matrix) , S- 
denotes a pseudo inverse matrix, and SI and S2 (column vectors) 
are echoes adjacent to the reference echo R (row vector) . 

Thereafter, the estimation coefficient A is used to 
estimate an unmeasured echo according to the expression (2) 
25 below (step 403) . 



12 



Su=AS' (2) 
where Su denotes an unmeasured echo, and S' equals [S'1,S'2] 
"T where S' 1 and S'2 denote echoes adjacent to Su. 

By performing the above processing, each echo is equally 
5 divided into Np parts as indicated with dashed lines 303 in 
Fig. 5. Fig. 5 is concerned with a case where Np equals 3. 
When an echo is divided, a more highly precise estimate is 
obtained than the one obtained when one echo is used as it is. 
An optimal number of parts into which one echo is divided is 

10 normally about seven . Assuming that the number of echo sampling 
points is one hundred and twenty-nine, one echo is divided into 
parts of, for example, eighteen, eighteen, eighteen, twenty-one, 
eighteen, eighteen, and eighteen sampling points. 

Moreover, the number of reference echoes should be equal 

15 to or larger than eight. For example, assuming that the number 
of reference echoes is four, the k-space is quartered. An angle 
6 occupied by one quarter is calculated as 360°/4 = 90°. In this 
case, echoes contained in the wide range of 90° have to be 
estimated using one estimation coefficient. An estimate of 

20 satisfactory precision cannot therefore be obtained. In 

contrast, when the number of reference echoes is eight or more, 
the range covering echoes to be estimated using one estimation 

coefficient is equal to or smaller than 45°. An estimate of 
satisfactory precision can therefore be obtained. The larger 
25 the number of reference echoes is, the higher the precision 
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in estimation is. However, a measuring time increases 
accordingly. Therefore, the number of reference echoes is 
normally set to about eight in terms of efficiency. 

Finally, measured echoes and unmeasured echoes that are 
estimated as mentioned above are combined and gridded. 
Thereafter, inverse Fourier transform is performed in order 
to reconstruct an image (step 404). 

Fig. 7 shows the results of the foregoing processing 
actually performed. Sixty-four echoes and eight reference 
echoes were measured, and fifty-six unmeasured echoes were 
estimated. The number of echo sampling points is one hundred 
and twenty-nine. Each echo is divided into seven parts, and 
estimation is performed. 

Fig. 7A shows an image resulting from the processing 
employed in the present invention. Fig. 7B shows an image 
resulting from reconstruction performed using sixty-four 
measured echoes alone. Fig. 7C shows an image resulting from 
reconstruction performed using one hundred and twenty-eight 
echoes. In Fig. 7B, numerous streak artifacts are seen in the 
background because of a small number of echoes. In contrast, 
in Fig. 7A, streak artifacts are hardly seen, and image quality 
is greatly improved. Consequently, the image quality is nearly 
equivalent to the image quality shown in Fig. 7C. 

As mentioned above, according to the present embodiment, 
part of unmeasured echoes is measured as a reference. An 
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estimation coefficient is calculated using echoes adjoining 
the reference, and the unmeasured echoes are estimated using 
the estimation coefficient. Consequently, since an echo 
measured as the reference is one of the unmeasured echoes, an 
5 imaging time hardly increases. Moreover, since the reference 
is used to calculate the estimation coefficient, unmeasured 
echoes can be highly precisely estimated compared when they 
are estimated by simple interpolation without using the 
reference . This results in a magnetic resonance imager capable 
10 of suppressing artifacts with little extension of the imaging 
time . 

( Second Embodiment ) 

In the pulse sequence shown in Fig. 1(A), a procedure 
starting with application of a slicing magnetic field gradient 
15 to measurement of an echo is repeated, for example, with a 
repetition time TR set to 4 ms, and the amplitudes of magnetic 
field gradients are changed so that one image can be produced 
using sixty- four echoes. 

In the present embodiment, the amplitudes of magnetic 
20 field gradients are changed so that the k-space will be scanned 

in a 9 direction at intervals of three echoes. Fig. 8(A) 
illustratively shows the scanning sequence. At this time, a 
frame rate is about 4 fps. A sliding window is adopted, and 
eight echoes are updated at a time. Consequently, the frame 
25 rate is increased to 32 fps. Thereafter, gridding is performed 
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in order to re-map echoes to a k-space 802 of a Cartesian grid. 
Finally, two-dimensional inverse Fourier transform is 
performed in order to reconstruct an image. 

Since the k-space 801 is scanned at intervals of three 
5 echoes, the frequencies of streak artifacts shift to higher 
frequencies. When the frequencies of streak artifacts shift 
to higher frequencies, the streak artifacts can be suppressed 
using a temporal filter (low-pass filter) . 

Furthermore, in order to suppress artifacts occurring 

10 on an edge of an image due to under-sampling, data to which 
the sliding window has not been applied and which is acquired 
by displacing scanning lines for fear scanning lines for 
adjoining frames may be superimposed on each other is employed. 
The temporal filter (low-pass filter) is applied to the 

15 time-sequential data, whereby artifacts derived from 
under-sampling can be suppressed. 

Fig . 9 shows the results 901 of decomposing a motion picture, 
which is reconstructed by scanning the k-space at intervals 
of three echoes, into frequencies in a temporal direction as 

20 an example of filtering of artifacts, and a temporal filter 
902 employed. In Fig. 9, the axis of abscissas indicates 
frequencies and the axis of ordinates indicates amplitudes. 
Reference numeral 903 denotes the frequency components of streak 
artifacts appearing in the motion picture . The temporal filter 

25 902 is applied to the frequency components 901 of the motion 
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picture, whereby the frequency components 903 of streak 
artifacts can be nullified. Consequently, the streak artifacts 
can be removed. 

Fig. 10 shows the results of imaging performed according 
5 to the present invention (first frames of motion pictures) . 
Fig. 10 (A) shows an image produced without scanning the k-space 
at intervals of three echoes, applying the temporal filter 
(low-pass filter) , and suppressing artifacts derived from 
under-sampling (scanning the k-space at intervals of no echo) . 

10 Fig. 10(B) shows an image produced by scanning the k-space at 
intervals of three echoes and applying the temporal filter 
(low-pass filter) . Streak artifacts are seen being suppressed. 
Fig. 10(C) shows an image produced by scanning the k-space at 
intervals of three echoes, applying the temporal filter 

15 (low-pass filter) , and suppressing artifacts derived from 
under-sampling. Not only streak artifacts attributable to a 
gap between adjoining data items but also artifacts appearing 
on an edge of an image are seen being suppressed. 

The present embodiment has been described on the 

20 assumption that the k-space is scanned at intervals of three 
echoes. At interval of how many echoes the k-space should be 
scanned is determined using the expression (3) below. 

(N+l)=Ne/2 (3) 
where N denotes the number of echoes N at intervals of which 

25 the k-space is scanned, and Ne denotes a magnification by which 
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a frame rate is increased by sharing echoes. 

The expression (3) is drawn out as described below. 
Assuming that a frame rate is multiplied by Ne by sharing echoes, 
a gap between adjoining data items varies in cycles of Ne/F 
5 (where F denotes a frame rate) . Accordingly, a streak artifact 
having a frequency F/Ne occurs. When the k-space is scanned 
at intervals of N echoes, the number of echoes to be scanned 
while a scanning line makes one turn is decreased to 1/(N+1) . 
Consequently, the frequency of a change in the gap between 
10 adjoining data items is multiplied by (N+1) . Accordingly, the 

frequency of a streak artifact comes to F/Nex(N+l). The 
expression (3) is drawn out under the condition of 

(F/Nex (N+1) =F/2) that the frequency of a streak artifact should 
be high (Nyquist rate F/2) . 

15 Fig. 11 describes an imaging procedure employed in the 

present invention. To begin with, parameters for an imaging 
sequence are entered in preparation for imaging (step 1101) . 

Thereafter, the N value representing the number of echoes 
N at intervals of which the k-space is scanned is determined 

20 according to the expression (3) (step 1102) . 

Thereafter, the sequence shown in Fig. 1(A) is applied 
in order to acquire echoes (step 1103) . At this time, as shown 
in Fig. 8, the k-space is scanned at intervals of N echoes (in 
the example shown in Fig. 8, N equals 3) . Moreover, when 

25 artifacts derived from under-sampling are suppressed, the same 



18 



number of echoes as the required number of frames is alternately 
scanned for fear scanning lines for adjoining frames may be 
superimposed on each other . 

Thereafter, when the number of echoes required for 
5 updating an image has been scanned, image reconstruction is 
performed (step 1104). 

Finally, the temporal filter (low-pass filter ) is applied 
in order to produce an image having streak artifacts suppressed 
(step 1105) . Moreover, when artifacts derived from 
10 under-sampling should be suppressed, the temporal filter 

(low-pass filter) is applied in order to suppress the artifacts . 

According to the present embodiment, the frequencies of 
streak artifacts derived from a gap between adj oining data items 
acquired by performing radial scanning is manipulated and 
15 filtered. Consequently, a marked advantage that the streak 
artifacts can be suppressed can be expected. 

Industrial Applicability 

According to the present invention, there is provided 
20 a magnetic resonance imager capable of efficiently suppressing 
artifacts in radial scanning. Moreover, a magnetic resonance 
imaging technology can be applied to an examination system or 
the like . This has significant meanings especially in the field 
of medicine. 

25 



